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AHSTRACT 


Rayleigh  waves  are  extracted  and  amplitudes  (with  standrrd 
errors)  and  origin  time  (within  a Cew  seconds)  determined  by  use  of 
reference  earthquakes  or  explosions,  high  pass  (HP)  and  normalized 
matcheo  (NM)  (or  correlation)  filters,  and  the  oceanic  group-velocity 
dispersion  between  38  and  18  seconds.  For  a signal-to-noise  ratio 
between  0.2  and  0.1  on  the  original  s ingle- component  records,  ampli- 
tudes for  Rayleigh  waves  over  oceanic,  paths  of  155°  at  station  MAT  and 
98°  at  station  KIP  have  been  determined  as  12  mg.  (10  digital  units) 
and  24  mu  (25  digital  units)  peak-to-peak  (p-p) , respectively,  with  a 
standard  error  of  less  than  117..  In  each  case  the  processed  signal  is 
the  highest  in  a half-hour  processed  record. 

The  method  also  provides  high  resolution  of  co-located  events  with 
short  time  separation,  or  widely  spaced  events  with  Rayleigh  waves 
arriving  nearly  simultaneously  at  a single  station.  Examples  include: 

(1)  clear  separation,  amplitude  determination,  and  so  forth,  at 
stations  KIP  and  MAT  of  two  Mg  =6.5  earthquakes  located  0.7°  and 
145  sec  apart  (off  the  coast  of  central  Chile);  (2)  clear  separation 
at  station  KIP  of  a Novaya  Zemlya  event  M,  = 4.8  from  interfering 
Rayleigh  waves  of  an  M.  = 5.0  Kermadec  Island  earthquake  arriving  120 
to  140  sec  prior  to  the  searched  event,  with  almost  complete  elimination 
of  interference  on  the  summed  vertical  and  radial  processed  components; 
and  (3)  clear  separation  at  station  KIP  of  two  co-located  Mb  = 4.4  and 
4.5  earthquakes,  6 minutes  apart,  from  the  area  off  the  coast  of  Chile, 
determination  of  their  amplitudes,  and  so  forth,  in  the  presence  of 
interfering  Rayleigh  waves  from  two  Central  Alaskan  earthquakes,  the 
first  (M  = 4.1)  arriving  15  minutes  prior  to  the  first  Chilean  Rayleigh 
wave  and^the  second  between  the  two  Chilean  arrivals. 


The  threshold  reached  (10  and  25  digital  units,  p-p)  for  stations 
MAT  and  KIP  at  98°  and  155°,  respectively,  corresponds  to  an  Mg  - 3.3 
and  probably  can  be  Improved  further. 
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INTRODUCTION 


Relative  amplitude  of  surface  waves  tliat  are  generated  for  a given 
short-period  P-wave  signal  (Press  et_  a_l_.  , 1963;  Brune  e_t  aK  , 1963; 
Ericsson,  1970;  Basham  and  Whitham,  1971)  is  one  of  the  most  useful 
discriminants  to  distinguish  earthquakes  from  explosions.  Short-period 
P waves  are  readily  detected  for  small  events,  but  for  small  events  the 
determination  of  the  surface  wave  magnitude  using  the  Rayleigh  waves  is 
much  more  difficult,  especially  when  the  signal  amplitudes  are  smaller 
than  the  noise  level  at  a single  ot  even  at  an  array  of  long-period 
stations.  Additional  difficulties  arise  when  interfering  surface 
waves  must  be  separated,  generated  either  from  multiple  arrivals 
(lateral  refractions  and  reflections  from  continental  margins),  from 
the  same  event,  or  from  other  events  at  nearly  the  same  location  and  a 
few  minutes  apart  or  at  distant  locations. 

Alexander  and  Rahenstine  (1967a,  1967b)  used  matched-filter 
techniques  to  improve  signal-to-noise  (S/N)  r- t io  for  LASA,*  LRSM,** 
and  observatory  stations  in  addition  to  pro  band-pass  filtering.  The 
method  compared  a known  Rayleigh  wave  signal  with  the  unknown  signal 
from  the  same  area  by  cross  correlation,  but  was  usually  applied  to 
seismograms  with  S/N  ratios  above  1 and  never  normalized.  The  beam- 
formed  matched-filter  S/N  ratio  for  13  LASA  stations  was  17  db  above 
the  mean  S/N  ratio  for  the  individual  band-pass  filtered  seismograms; 
for  13  LRSM  stations  it  was  15  to  16  db  above  the  mean  S/N  ratio  of 
band-pass  filtered  seismograms.  However,  Alexander  and  Rabenstine's 
approach  (1967b,  p.  14)  did  not  allow  reliable  signal  amplitude 
estimates  where  signals  were  below  the  highest  noise  peaks;  but  they 
showed  from  synthetic  records  that  an  S/N  ratio  of  0.35  still  allowed 
extraction  of  the  teleseismic  surface  waves.  This  report  presents  a 
reliable  method  of  determining  surface  wave  signal  amplitude  (and 
standard  error)  that  works  at  S/N  ratios  between  0.2  and  0.1)  for 
single  seismograms. 

McDonald  et  al^  (1974)  extended  the  Alexander-Rabenstine 
approach  by  synthesizing  a matched  filter  with  the  same  average 
group-velocity  dispersion  as  three  Greenland  Sea  and  one 
Chinese  events,  as  recorded  at  Grand  Saline,  Texas.  They  applied  the 
filter  successfully  to  the  small  earthquakes  in  the  Sino-Soviet  land 
mass,  taking  full  advantage  of  the  long-period  surface  wave  guide  over 
the  pole.  The  wave  guide  extends  to  periods  of  about  75  sec  and 
generates  a well-dispersed  wave  train.  A detection  threshold  of 
M = 3.5  was  obtained  for  the  Rayleigh  waves  arriving  at  the  single 
vlrtical  component  at  Grand  Saline  llnough  the  polar  wave  guide,  an 
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M,  = 4.5  equivalent  for  focal  depth  less  than  50  km.  Their  use  of  the 
chirp  filter  between  75  and  20  sec  made  the  detection  in  the  wave 
guide  beam  (about  45°)  independent  of  a specific  reference  event,  but 
did  not  seem  to  allow  separation  of  multiple  events  from  multiple 
arrivals  (from  a single  event),  nor  *-o  provide  amplitude  information 
from  the  correlation  processing. 

Capon  et_  al.  (1969)  and  Capon  (1970)  successfully  increased  the 
S/N  ratio  for  Rayleigh  waves  (at  1ASA)  using  noise  prediction,  matched 
(chirp)  filtering,  and  beam-forming  in  the  frequency -wave  number  domain. 
Later,  Capon  and  Evernden  (1971)  applied  this  method  to  the  40-sec 
Rayleigh  waves.  Only  the  vertical  components  were  used,  since  the 
horizontal  instruments  were  too  noisy  to  apply  the  surface  wave 
discriminant  at  low  S/N  ratios.  Capon  et^  aJL..  (1969,  table  2)  obtained 
total  average  S/N  ratio  gains  of  22  db  (for  matched  chirp  filter,  noise 
prediction  filter,  and  beam-forming  for  17  seismometers).  Subsequently, 
Capon  and  Evernden  (1971)  considered  the  resolving  limits  of  the  method 
(at  40  sec)  as  to  azimuth  and  time  separation,  and  found  that  events 
from  almost  opposite  directions  (such  as  from  the  Kuril  Islands  and 
from  Ecuador),  simultaneous  arrivals,  and  power  level  differences  of 
6 db  (average)  could  be  separated.  With  increasing  time  difference, 
larger  level  differences  could  bo  resolved  with  an  average  relative 
amplitude  of  between  26  and  30  db  at  a time  separation  of  more  than 
10  minutes,  the  actual  maxima  and  minima  of  the  resolvable  level 
differences  depending  on  the  die-off  character  of  the  40-sec  Rayleigh 
wave  (Capon  and  Evernden,  1971,  Fig.  12). 

While  the  use  of  reference  earthquakes  (or  explosions)  requires 
more  detailed  information  than  does  the  chirp  matched-filter  method, 
the  advantages  of  the  former  are  that  the  signal  amplitude  (and  its 
standard  error)  can  be  obtained.  Multiple  arrivals  from  the  same 
source  with  different  travel  paths  are  included  in  the  (earthquake 
reference)  matched  filtr,  and  therefore  do  not  appear  as  distinct 
arrivals  on  the  correlation  trace.  The  matched  chirp  filter  however 
shows  these  multiple  arrivals  and  is  unable  to  distinguish  them  from 
multiple  events.  In  addition  use  of  a normalized  correlation  (=  1 for 
perfect  signal  match)  permits  the  specification  of  decision  levels, 
recognition  of  defective  recording  traces  and  other  erroneous 
signatures  that  are  not  readily  computerized,  and  application  of 
statistical  measures.  As  an  example,  the  Rayleigh  wave  for  an 
Mi,  = 4.5,  Ms  = 3.3  Chilean  earthquake  has  been  extracted  at  the  MAT 
(A  = 17000  km)  high  gain  long-period  vertical  component  (after 
preliminary  high-pass  filtering);  the  amplitude  was  determined  as 
9.4  digital  units  (peak- to-peak)  (about  12  my,  at  30  sec).  The 
correlation  peak  is  the  highest  in  the  30-minute  period  scanned.  The 
noise  background  at  the  tHe  was  between  90  and  100  digital  units  (p-p) 
on  the  original  unfiltered  record  (see  Outputs  4 and  4/40  at  end  of 
report) . 

In  this  study,  a system  tic  attempt  was  made  to  determine  the 
influence  of  the  Length  of  the  reference  record,  and  of  the  high-pass 
prefiltering. on  the  results  of  single  and  combined  vertical  and  horizon 
tal  component  processing.  The  best  detection  and  separation  from  a 


single  station  were  achieved  by  summing  the  vertical  and  radial  Rayleigh 
wave  components  after  40-sec  high-pass  and  mat -lied  filtering  (in  that 
order).  However,  so  far  we  have  considered  mostly  oceanic  paths,  and 
have  taken  advantage  of  the  well-dispersed  oceanic  Rayleigh  waves  for 
periods  shorter  than  40  see,  the  worldwide  noise  minimum,  and  the 
maximum  gain  of  the  high  gain  long-period  stations  between  30  and 
40  sec  (Murphy  et  a 1 . , 1972). 

NORMALIZED  MATCHED  FILTER 

The  source,  transmission  path,  and  instrument  response  function 
determine  the  seismic  signal  recorded  on  a seismograph.  When  using  a 
reference  earthquake  as  the  matched-filter  impulse  response,  it  is 
tacitly  assumed  that  the  source  function  (amplitude  and  phase)  can  be 
simply  scaled  up  or  down  (the  focal  mechanism  remaining  more  or  less 
constant)  while  the  path  and  instrument  functions  remain  constant.  Tni s 
differential  approach,  which  eliminates  the  need  for  detailed  knowledge 
of  the  transfer  functions  of  the  transmission  path  and  receiver,  has 
been  applied  during  other  irvest igat ions  (Berckhemer,  1962;  Alexander 
and  Rabenstine,  1967a  and  b;  and  lately  by  Weidner  and  Aki,  1973). 

Several  papers  (Berckhemer,  1962;  Aki,  1967;  Savage,  1Q72;  Haskell, 

1964;  Bruno,  1968)  give  theoretical  fault  models  and  experimental  data 
for  the  spectral  amplitude  density  as  a function  of  fault  parameters, 
such  as  rupture  pattern  in  the  fault  surface,  elastic  and  rupture 
velocities,  and  fault  dimensions.  All  these  investigations  show  that 
for  3 < Ms  < 5.5,  the  source  spectrum  is  constant  for  periods  longer 
than  20  sec  (Aki,  1967,  Figs.  3 and  4;  or  Berckhemer,  1962,  Fig.  lid). 

We  therefore  feel  justified  in  assuming  a constant  spectrum  requiring 
. only  an  amplitude  scaling  of  the  Rayleigh  waves  to  be  extracted  for  the 

earthquakes  (of  similar  focal  mechanism,  area  and  depth),  as  long  as 
the  reference  earthquake  does  not  exceed  Ms  = 5.5.  On  the  other  hand, 

M =5.5  assures  a healthy  S/N  ratio  for  the  reference  earthquake  on  the 
high  gain  long-period  station  records  if  the  Rayleigh  wave  is  not 
contaminated  by  other  surface  wave  trains.  The  assumption  of  similar 
focal  mechanism  is  borne  out  by  the  fact  that  the  correlation  peaks 
always  show  a positive  sign. 

The  basic  formulas  for  the  correlation  traces  are  given  in  the 
Appendix;  they  present  the  cross  correlation  coefficient  of  the  scanned 
record  trace  with  respect  to  the  reference  trace  (Crow  e£  al . , I960). 

A matched-filter  response  (after  high  pass  filtering  at  40  sec)  is 
presented  in  Figure  2b  for  the  Z component  of  the  reference  earthquake 
in  the  Kurils  (Output  2/40).  A perfect  match  would  result  in  an  output 
of  "one"  at  the  lag  time  of  the  beginning  of  the  reference  signal.  The 
advantage  of  the  normalized  matched-filter  output  is  that  it  permits 
determination  of  an  absolute  value  for  the  maximum  noise  correlation 
with  the  reference  signal  (of  a given  length  and  area)  which  value,  if 
exceeded,  indicates  the  presence  of  a seismic  (surface)  wave.  The 
correlation  is  formed  between  the  reference  and  the  record  by  increment- 
ing lags.  When  a maxi^'m  correlation  is  found,  this  correlation  time 
is  taken  as  the  start  l .me  of  the  hidden  signal  in  the  record  (as 
plotted  in  the  correlation  trace  above  the  record  trace).  To  make  the 
maximum  correlation  sta  on  independent  for  a given  earthquake-station 
combination,  the  time  labeling  (on  the  correlation  trace)  is  shifted 


Fig.  1 


One-third  the  sum  of  vertical,  N-S,  and  E-W  correlation  t 
function  of  data  length  for  the  Kuril  Islands  earthquakes 
recorded  at  KIP.  For  99  per  cent  level,  refer  to  te  . 
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(to  an  earlier  time)  by  the  difference  between  the  origin  time  and  the 
start  time  of  the  reference  earthquake  trace  . This  insures  that 
irrespective  of  the  station  used,  the  maxirium  correlation  occurs  near 
the  scanned  earthquake's  origin  time,  assuring  proper  identification 
at  several  stations.  Details  of  all  earthquakes  and  explosions  processed 
or  mentioned  are  pres  *nted  in  the  table  or:  page  19. 

An  example  from  station  MAT  is  given  in  Output  1.  The  fourth, 
seventh,  and  last  traces  are  the  Rayleigh  wave  portions  of  the  Z,  N-S, 
and  E-W  components  of  the  reference  earthquake.  The  third,  sixth,  and 
ninth  traces  are  the  scanned  Z,  N-S,  and  E-W  records  (note  that  the 
E-W  component  was  not  operating  properly,  since  no  signal  is  visible). 
The  second,  fifth,  and  eighth  traces  are  the  normalized  matched  filter 
outputs  with  scale  from  +1  to  -l,  and  the  "SM"  (first)  trace  is  one- 
third  the  sum  of  the  three  matched-filter  traces.  If  only  two 
components  have  been  processed  ('  'ch  as  vertical  and  radial)  the  SM 
trace  is  one-half  the  sum  of  the  two  correlation  traces.  Figure  1 
shows  the  signal  correlation  coefficient  obtained  (for  the  same  Kuril 
earthquakes  recorded  at  KIP)  as  a function  of  the  drta  lengths  of  the 
reference  Rayleigh  wave  train,  and  the  three  or  four  highest  (positive) 
noise  peaks  (for  a half-hour  scan).  It  is  clear  that  they  follow  the 
statistical  distribution.  The  99  per  cent  line  means  that  only  1 per 
cent  of  the  population  of  correlation  coefficients  with  values  higher 
than  indicated  by  that  line  are  associated  with  data  sets  (six 
variables:  three  reference  traces  and  their  noise  and  three  scanned 

traces  and  their  noise)  that  do  not  show  any  relation.  On  the  other 
hand,  the  decision  level  (to  consider  a correlation  level  as  presenting 
a signal)  depends  on  the  reference  data  length.  After  this  preliminary 
investigation,  w:  always  used  an  841  data  point  (840  sec)  reference 
section.  These  14  minutes  of  reference  section  correspond  to  a variable 
portion  of  the  group-velocity  dispersion  curve  depending  on  station  to 
epicenter  distance,  but  seem  to  yield  good  results  for  extraction  and 
S/N  ratio  of  the  correlation,  since  for  the  oceanic  path  they  include 
the  most  dispersed  portion  of  that  curve  up  to  over  10,000  km  distance. 

Since  the  relatively  long  reference  section  includes  at  least  the 
early  multiple  arrivals  (if  they  exist),  they  all  will  be  condensed 
into  the  correlation  peak  to  help  separate  it  from  other  events.  This 
is  a distinct  advantage  of  the  reference  earthquake  matched-filter  over 
the  chirp-type  filter  (that  usually  shows  the  multiple  arrivals  as 
separate  signals)  in  extracting  small  surface-wave  events  (as  well  as 
amplitude  determinations) . 
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AMPLITUDE  DETERMINATION 


Another  disadvantage  of  the  chirp-filter  is  its  complete  lack  of 
attenuation  and  instrument  response  information,  so  that  amplitude 
scaling  is  unobtainable.  The  reference  earthquake  matched-filter,  by 
contrast,  contains  all  this  information  and  therefore  is  ideally 
suited  for  amplitude  extraction.  In  addition,  all  the  computation 
steps  required  to  determine  the  correlation  coefficient  between  the 
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Fig.  2B.  Power  spectrum  of  the  matched-filter  time  function  of  the  KIP-Z 

Kuril  Islands  reference  earthquake  of  Output  2/40  (not  normalized) 
Note  the  reference  also  was  HP,  prefiltered  at  40  sec. 


scanned  and  reference  trace  contain  all  the  information  for  determination 
of  the  relative  amplitude  and  the  standard  error,  since  a simple  linear 
relation  between  reference  earthquake  and  scanned  earthquake  is  assumed. 
Determining  the  relative  amplitude  involves  finding  the  slope  of  the 
regression  line  of  tne  scanned  data  (at  time  of  maximum  correlation)  on 
the  reference  data  and  the  standard  error  of  the  slope.  (See  Crow 
et  aj_.  , I960,  p.  152) . 

Since  the  relative  standard  error  of  the  slope  is  only  a function 
of  the  reference  length  used  and  the  correlation  obtained,  the  decision 
on  an  extraction  level  the  presence  of  a signal  at  a predetermined 
value  of  correlation)  also  sets  the  upper  limit  for  the  standard  error 
in  per  cent  of  the  amplitude.  (Figure  7 gives  the  relative  standard 
error  as  a function  of  the  correlation  coefficient  obtained  for  a 
reference  length  of  841  data  points;  the  corresponding  formulas  are 
presented  in  the  Appendix.)  Tli*s  information  (time  of  maximum 
correlation,  slope  * relative  amplitude,  standard  error)  was  only 
printed  out  for  the  maximum  positive  correlation  in  each  (usually  half- 
hour)  scan  and  was  added  at  the  end  of  the  corresponding  correlation  in 
the  output  plots. 

Since  for  nuclear  explosions  and  their  cavity  collapse  (and  possibly 
aftershocks)  a negative  correlation  is  expected  (and  observed,  see 
Output  8R)  , amplitude  and  error  information  should  be  obtained  for  the 
minimum  correlation  as  well. 

Simply  taking  the  log10  of  the  slope  (or  the  relative  amplitude) 
gives  the  difference  in  Ms  between  the  extracted  and  the  reference 
earthquake.  Since  Ms  for  the  reference  earthquake  can  be  determined  as 
an  average  over  many  stations  (or  may  include  considerations  for 
radiation  patterns),  the  value  determined  for  the  scanned  extracted 
earthquake  is  free  of  individual  station  corrections  and  appears  to  be 
very  stable  (see  Fig.  6). 


PREFILTERING 

Most  reference  earthquakes  (over  the  oceanic  path  used  so  far) 
showed  the  longest  periods  between  35  and  39  sec,  and  the  noise 
appears  large  at  longer  periods,  especially  on  the  horizontal 
components  of  the  high  gain  system.  Therefore,  high-pass  filters  with 
cutoff  periods  of  60,  50,  40,  and  30  sec  have  been  applied  before 
matched  1 iltering  to  determine  improvement  in  correlation  and  the 
S/N  ratio  (filter  characteristics  are  presented  in  Fig.  2A  and  2B) . 

The  value  used  here  is  the  (correlation  signal  peak) /(next  highest" 
positive  noise  peak)  (C:SN)  for  a hull -hour  scan  and  a reference  length 
of  840  sec,  correlated  after  both  the  reference  and  the  scanned  trace 
have  been  high-pass  filtered. 

The  general  relation  between  the  correlation  and  C:SN  is 
depicted  in  Figure  3 for  the  earthquake  from  off  the  coast  of  Chile 
recorded  at  MAT  and  KIP  and  t h,  Kuril  Islands  earthquake  recorded  at 


CORRELATION 


Fig.  3.  Correlation  coefficient  (C)  versus  correlation  signal  noise  (C:SN) 
for  various  HP  filters.  Increase  in  correlation  results  in 
improved  amplitude  determination;  increase  in  C:SN,  in  improved 
detection.  Pre  high-pass  filtering  near  40  sec  to  50  sec  gives 
best  results  for  pure  oceanic  travel  path. 
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KIP  (see  table,  p . 18).  Underlined  data  points  are  values  obtained  without 
filter.  Hie  progress  along  a line  corresponds  to  IIP  filters  at  b(),  50, 

AO,  and  30  sec,  except  for  the  MAT-Cliile  Z trace  where  the  last  point 
corresponds  to  a AO-sec  HP  filter.  It  should  be  noted  that  placing  a 
decision  level  for  the  presence  of  a signal  at  0.3  for  the  correlation 
of  a single  trace  (with  8A0-sec  reference  in  ; half-hour  scan)  seems  to 
assume  that  the  highest  signal  peak  is  not  due  to  noise  (or  C:SN  > 1 for 
all  signals  with  C > 0.3).  Note  also  how  the  progressive  filtering 
pulls  the  earthquake  signal  out  of  the  noise  for  the  Chilean  earthquake 
recorded  on  KIP  Z,  E-W,  and  MAT  Z.  Outputs  2 through  A/AO  show  most  of 
the  corresponding  records.  The  MAT  Z-amplitude  of  the  Rayleigh  wave 
(P*P)  after  the  AO-sec  HP  filtering  was  9. A digital  units  (12  ni|i,  p-p) , 
whereas  the  noise  on  the  original  unfiltered  Z was  about  90  to  100 
digital  units  (p-p)  (see  Output  A).  Figure  A presents  the  same  single 
component  data  as  Figure  3 plus  those  for  the  rotated  radial  horizontal 
component  (R) , and  the  correlation  stun  traces  of  either  the  three 
components  (SM)  or  the  radial  and  vertical  components  (RSM)  tor  the 
optimum  HP  filter.  It  is  clear  (Fig.  AA)  that  correlation  increases  as 
more  and  more  of  the  long-period  noise  is  filtered  out  (to  A0  sec) 
until  part  of  the  signal  becomes  filtered  out  as  well,  it  least  for 
weak  signals  such  as  those  from  the  Chile  earthquake.  The  same  observa- 
tion holds  for  the  C:SN  (Fig.  AH).  Increased  correlation  results  in  a 
smaller  standard  error  tor  the  amplitude  determination,  increased  C:SN, 
in  better  detection.  It  is  also  seen  that  the  radial  component  alone 
gives  both  higher  correlation  and  better  C:SN  than  either  of  the 
horizontal  components.  In  effect,  this  is  the  same  as  polar izat ion 
filtering.  The  same  improvement  also  holds  for  the  rotated  sum  (RSM) 
trace  over  the  SM  trace  for  KIP.  It  was  not  possible  to  obtain 
rotated  data  for  MAT  since  the  E-W  component  was  not  operating  properly. 

Most  of  the  individual  output  plots  from  which  the  data  points  were 
obtained  are  presented  in  the  output  section.  It  should  be  stressed 
that  all  Rayleigh  wave  signals  were  well  In  low  the  noise  level  of  Un- 
original records  (Outputs  2,  3,  and  A). 


If  S/N  ratio  of  the  original  record  is  deliind  as  (p-p)  filtered 
s igna I ) / (p-p  unliltered  noise),  all  in  digital  units,  then  C:SN  is 
the  S/N  ratio  ol  the  processed  record;  an  improvement  would  he  (C:SN)/(S/N 
original  letord).  Improvement  ol  the  normalized  matched  filter  for  KIP 
without  lll'-t  i I tering  on  a single  trace  seems  to  be  near  12  db  .md  at 
about  lb  to  19  db  when  AO-sec  HP-liltering  is  applied.  However,  very 
weak  signals  (as  Z MAT- -Chile  A0-sec  IIP)  seem  to  show  higher  improvement 
(22  db).  This  Is  partially  due  to  tin  limit  on  tin-  eorre  I .it  i on  i,. un- 
being = I (or  a C : SN  of  nearly  5 sine  tmise  correlated  with  th 
reference  signal  hy  about  0.2  or  le-,.;-. 
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AMPLITUDES  AND  RELATIVE  MAGNITUDES 


Figure  5 presents  the  umpl  iti.jo  (p-p)  determinations  from  the 
variously  filtered  outputs.  The  individual  amplitudes  are  obtained  by 
multiplying  the  p-p  amplitude  of  the  reference  earthquake  by  the  relative 
amplitude  (slope)  of  the  scanned  earthquake.  As  an  example,  consider 
the  vertical  component  of  the  Kuril  Islands  earthquake  of  October  3 
(day  276),  1973  at  station  MAT  (in  Output  1).  Maximum  correlation 
occurs  where  the  start  of  the  reference  signal  is  positioned  at  the 
vertical  line  (near  13H01M).  The  vertical  correlation  trace  time  label 
is  displaced  to  an  earlier  time  by  the  difference  between  the  origin 
time  and  start  time  of  the  reference  earthquake  Rayleigh-wave  trace. 

The  peak  occurs  at  12H54M52Scc  (pr  inted  out  as:  "i  me  of  maximum  is: 

1973-276-12:54:52"  on  the  right-hand  side).  This  time  is  the  origin 
time  (see  table,  p.  19).  The  reference  earthquake  had  a maximum  p-p 
amplitude  of  4585  digital  units  (from  the  lower  edge  to  the  upper  edge 
of  the  left-hand  scale,  repeated  over  the  word  "station  for  clarity). 
This  value  multiplied  by  the  relative  amplitude  of  the  signal  in  the 
scanned  trace  gives  the  maximum  p-p  signal  amplitude  of  the  correlated 
earthquake.  The  relative  amplitude  is  printed  out  on  the  right-hand 
side  under:  "Slop,  is  0.069446".  Therefore,  the  amplitude  of  the 

vertical  signal  for  this  Kuril  Islands  earthquake  recorded  at  MAT  is 
318.4  digital  units  (with  a standard  error  of  ± 4585  • 0.0016987  or 
± 7.8  digital  units).  Note  that  usually  the  maximum  p-p  amplitude  on 
the  record  (356  digital  units  near  13H04M30Scc)  would  be  picked  to 
determine  the  magnitude  from  this  record.  The  log10  of  the  relative 
amplitude  (logl0  0.069446  « - 1.16)  gives  the  difference  in  surface- 
wave  magnitude  (-  1.16)  between  the  Kuril  Islands  extracted  and 
reference  earthquakes.  Similarly,  the  difference  in  magnitude  at 
station  KIP  (see  vertical,  Output  2)  is  log^Q  0.071356  = - 1.15. 

Figure  6 indicates  that  the  magnitude  difference  AMS  is  a stable  value 
at  a given  station,  and  that  even  for  the  extremely  weak  signals 
extracted  for  the  Chilean  earthquake  the  variation  is  only  about  0.1. 
Since  the  relative  standard  error  SE(b)/b  ol  the  slope  b (b  - relative 
amplitude  with  respect  to  the  reference  earth ;uake)  is  nothing  more 
than  a function  of  reference  data  length  and  the  correlation  obtained 
(see  Appendix  for  derivation),  a single  curve,  is  valid  for  all 
correlation  outputs  with  a fixed  reference  length.  Figure  7 gives  the 
standard  error  in  per  cent  of  the  amplitude  as  a function  of  the 
correlation  and  for  N * 84!  data  points  (or  840  sec). 


The  trend  in  Figure  5 indicates  that  one  can  expect  to  see  signals 
with  about  3 digital  units  p-p  a a correlation  level  of  0.2  or  so,  a 
signal  that  could  be  useful  if  properly  identified.  The  identification 
can  have  several  forms.  Such  a s'gnal  level  still  would  be  determined 
with  a standard  error  in  amplitude  of  less  than  20  per  c nt  (Fig. 

With  such  low  amplitudes  of  a few  di|ital  units,  the  limits  of  the 
recording  system  might  be  reached,  rather  than  those  of  the  method. 
Sine'  the  oceanic  path  signals  from  the  Kuril  Islands  and  the  area  off 
the  Chile  coast  show  the  maximum  p-p  record  amplitude  near  30  sec, 

3 dig  til  units  (r  3 
Mr  (30)  - 2.4  at  100° 


1/3  mu)  will  place  the  detection  threshold  noa- 
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HIGH  PASS  FILTER  PERIOD 


Fig.  4a. 


Fig.  4b. 


Correlation  as  a function  of  high  pass  filter  period.  Improved 
correlation  results  in  better  signal  amplitude  determination. 

SM  is  one-third  the  sum  of  the  vertical  N-S  and  E-W  component 
correlation  and  RSM  is  one -ha If  the  sum  of  the  vertical  and  radial 
component  correlation.  R is  the  radial  component  correlation. 

Note  that  RSM  is  a considerable  improvement  over  SM  and  similarly 
the  correlation  of  the  radial  component  is  better  than  either  E-W 
or  N-S  components.  Note  that  C = 0.25  seems  to  be  the  correlation 
where  the  signal  emerges  as  the  highest  correlation  peak  in  the 
half-hour  of  scanning. 

Earthquake  correlation  peak/highest  positive  correlation  noise  peak 
on  the  half-hour  record  as  function  of  HP  filter  period.  For 
C:SN  > 1 the  earthquake  signal  is  the  highest  on  the  half-hour 
record.  The  maximum  in  the  s igr.al/noise  corresponds  to  the  maximum 
gain  of  the  Instruments  and  the  suppression  of  noise  in  the  period 
range  bevond  40  sec,  where  no  signal  is  observed  for  the  pure 
ocoanic  path. 
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SIGNAL  AMPLITUDE  o MAT 

{ P-P ) DIG  UNITS  - KIP 


Fig.  5.  Signal  amplitudes  extracted  by  normalized  matched  filtering  after 
high  pass  filtering,  at  periods  indicated,  for  single  component 
records.  In  all  cases  the  correlation  (positive)  peak  was  the 
highest  during  the  half-hour  scan.  Standard  error  of  the  amplitudes 
with  a correlation  (C)  near  0.3  is  about  11  per  cent  and  decreases 
to  6 per  cent  for  C near  0.5. 
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HP  FILTER  PERIOD 


Di  f ference  in  surface  wave  magnitude-  of  extracted  signal  with 
respect  to  reference  earthquake . Note  the  very  small  variation  for 
various  HP  filters  or  between  MAT  and  KIP  for  th  unfiltered 
records  (Kuril  Islands  earthquake).  In  all  records,  signal  level 
was  below  noise  (except  MAT-Kur i 1 ) ; in  the  case  of  MaT-Z  Chile,  by 
as  much  as  a factor  near  10. 
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SEPARATION  OF  OVERLAPPING  WAVE  TRAINS 


iAt  a single  si  a'  ion  the  (earthquake  reference')  matched  filter  also 
separates  out  either  co-located  events  with  small  time  differences,  or  a 
widely  separated  (undesired)  event  when  the  Rayleigh  wave  arrivals  are 
almost  simultaneous,  or  both.  The  relative  amplitudes  can  be  obtained 
as  exp  1 a i ned  ear  1 i or . 

Outputs  5,  6,  and  7 show  the  separation  of  two  M ■ 6.5  earthquakes 
off  the  coast  of  Chile  that  occurred  l*+5  sec  and  about  0.7°  apart.  A 
very  clear  separation  is  obtained  at  KIP  (Output  5)  and  MAT  (Output  6) 
despite  the  relatively  noise-contaminated  reference  used  for  MAT.  Even 
the  MAT  SM  correlation  trace  shows  the  separation  despite  the  one-third 
contribution  from  the  non-operating  E-W  component.  Note  also  that 
neither  rotation  nor  filtering  was  applied.  Output  7 is  tht  ore-half 
(SM  KIP  + SM  MAT)  s ignal,  c lear ly  defining  the  signals  as  originating  in 
the  area  off  the  coast  of  Central  u^ile. 

Output  8R  shows  an  example  of  extraction,  at  KIP,  of  a Novaya  Zemlya 
event  Mb  = 4.8  with  negative  correlation.  An  interfering  Rayleigh  wave 
arrived  about  2 minutes  earlier  (approximate  arrival  underlined)  from 
the  Kermadec  Islands  (Mb  = 5.0).  Since  an  explosion  was  used  as  a 
reference  trace,  the  negative  correlation  clearly  indicates  that  the 
later  event  was  not  an  explosion  but  rather  an  earthquake  or  a cavity 
collapse.  Note  that  the  interfering  Rayleigh  wave  (from  the  S-W  of 
KIP)  has  a dispersed  Z component  correlation  that  is  quite  different 
from  the  negative  correlated  signal,  and  that  the  SM  trace  nearly 
eliminates  this  undesirable  signal.  SM  presents  nearly  the  rotated  SM 
output,  since  the  signal  from  Novaya  Zemlya  arrived  approximately  from 
the  north. 

Finally,  Output  4 shows  the  original  data  for  two  earthquakes  from 
off  the  coast  of  central  Chile  (Mb  = 4.4  and  4.5)  , separated  by 
6 minutes.  Interfering  Rayleigh  waves  from  two  Central  Alaska  earthquakes 
arrived  15  minutes  before  the  first  one  and  just  between  the  two 
(approximate  arrival  time  underlined  on  Output  9/40R,  the  hp-filtered 
and  rotated  record  for  these  same  earthquakes).  The  computer  correctly 
picked  one  of  them  on  the  vertical  and  the  other  on  the  radial  trace. 
Signal  peaks  on  the  RSM  trace  are  clearly  the  highest  during  the 
processed  time  interval. 


E Earthquakes  (USGS  Monthly  Listings,  October  1973) 
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APPENDIX 

"0"  Phase-Shift  High  Pass  Filter 

If  A ( i ) is  the  original  data  trace  j = 1,2... 

digitized  at  a t ime  interval  At, 
the  filtered  output  AF(j)  is  given  by: 

2N+1 

AF( j)  = £ A ( j -1 -N+i)  ■ F(i) 
i-l 

and  F (i)  = T(i)  • W(i) 

where  T(i)  = tapering  function 

W(i)  = Fourier  transform  of  the  frequency  window. 
The  Hamming  tapering  function  is: 

T ( i)  = 0.54  + 0.46  cos  [ ] 

and  W(i)  for  a high  pass  filter  is: 

sin  [ (-N-l+i)n2AtfH  ] 


W(i)  = 6 (N+l-i)  - 2AtfH 


(-N-l+i)n2AtfH 
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Rotat i on 


The  rotation  program  is  applied  to  either  the  direct  data  or  the 

filtered  data.  Only  a fixed  average  gain  factor  is  used;  however,  at 

present  the  gain  factor  is  variable  with  frequency. 

A3  = amplitude  of  Z component  in  digital  units 
A2  = amplitude  of  N-S  component  in  digital  units 

A3  = amplitude  of  E-W  component  in  digital  units. 

If  cp  = angle  at  station  clockwise  from  North  to  epicenter,  the  trans- 
formation is  given  by: 


A l « A 1 

R = sin  cp  i-  A2  cos  cp  (radial) 

T - G3A3  (-  cos  cp)  t 1^2  sin  cp  (transverse). 

The  gain  factor  G3  is  set  so  that,  as  an  average  over  the  predominant 
frequency  range, 


A2 (digital  units/p  of  calibration 
^3  A3 (digital  units/p  of  calibration 


The  transformation  gives  the  output  R , T in  digital  units,  approximately 
at  the  gain  of  component  A2(N-S)  with  R positive  towards  the  epicenter. 
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Correlat ion 


The  cross  correlation  (normalized  matched  filtering)  between  a 
reference  earthquake  and  a scanned  record  in  performed  either  for 

(a)  the  direct  data  as  recorded  and  for  each  individual  channel 
(Z,  N-S,  E-W); 

(b)  the  high-pass  filtered  data  with  the  same  filter  applied  to 
the  reference  earthquake  and  the  scanned  data  first;  or 

(c)  the  rotated  (or  rotated  and  filtered)  data  and  for  the 
channels  Z,  AR,  AT. 

The  following  gives  the  cross  correlation  for  a single  component 
(as  under  a,  b,  or  c (At  = 1 sec  being  omitted  fran  the  formulas)): 


C(j)  - C(t)  - C(ts  - 1 + j - (tgg  “ 

tD„  = origin  time  of  reference  earthquake  (from  bulletins, 
KU 

and  so  forth) 


r 
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N N 


N^R(i)  S (tg-l+j+i)  -^R(i)  * J S(ts-Hj+i) 


C(t)  = 


i = l 


i=l  1-1 


[ N y R 2 ( t ) ~ ( ^ R(i) ) 2 ]'"  [ N ^ S2  (t  s- 1+j+i)  - ( ^ S(ts-l+j+i))2 
i=l  i=L  1=1  i =1 


j = l,2,...(tE-ts+l) 


/TT  / in 


Standard  Error 


The  slope  of  the  regression  line  of  the  scanned  data  S(j+i)  on  the 
reference  data  R ( 1 ) Is  given  by: 

b(t)  = slope  (j)  = , 


and  the  standard  error  of  b(t)  Is: 


SE(b) 


( 111 
( t r 


- b ( t) > 


2 


In  the  present  program,  the  computer  retains  and  prints  out  the  value  of 
b and  Its  standard  error  only  for  the  time  of  maximum  correlation.  The 
slope  Is  then  the  relative  amplitude  of  the  extracted  earthquake  with 
respect  to  the  reference  earthquake  If  the  former  was  correctly  picked. 
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Since  “ T 

C 


ir 


ii  • hi  hi 
i2  ii 


SE(b)  ■ [ sr2  (|z  - b2  ) ] 


or  the  relative  standard  error  of  b 


is  only  a function  of  reference  data  length  and  the  correlation 
obtained . 


■ 
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Explanation  of  Output  Plots 


Output  numbers  are  for  reference  in  text  and  table  of  earthquakes. 
The  same  number  is  used  for  the  same  station  - scanned  earthquake 
combination.  A number  (60,  50,  40,  or  30)  after  a virgule  (/)  Indicates 
HP-filtered  data  at  that  particular  period  for  both  the  reference 
earthquake  and  the  scanned  record;  R indicates  rotated  data  for  the 
horizontal  components. 

Trace  designation  Ap  A£,  original  or  HP  filtered  Z,  N-S,  E-W; 

Cp  C2>  C3:  correlation  of  Aj  with  Aj  reference,  and  so  forth; 

SM:  (sum  of  C traces)/ (number  of  correlation  traces),  R,  CR:  radial, 

radial  correlation. 

Start  time  is  given  after  station  identification.  Tick  marks  are 
2 minutes  apart.  Amplitudes --(on  the  left)  in  digital  units  (not 
readable);  p-p  in  digital  units,  was  redrawn  over  the  word  "STATION" 
in  most  records. 


Output  # 1,  2,  2R,  2/60,  2/50,  2/40,  2/30  Kuril  Islands 

3,  3/60,  3/50,  3/40,  3/40R,  3/30  1 Off  the  coast  of  Central 

4,  4/40  J Chile,  extraction 

5,6,  7 Off  the  coast  of  Centre! 

Chile,  separation 

Novaya  Zemlya, 
extraction,  separation 

Off  the  coast  of  Central 

Chile,  extraction  and 

separation. 

K.  * 4.5  and  4.4 
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Preceding  page  blank 
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